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Abstract. We review our recent results on exclusive production of fj,'^ , heavy quark- 
antiquark, and meson-antimeson pairs in ultraperipheral, ultrarelativistic heavy ion collisions. 



1. Introduction 

Ultrarelativistic collisions of heavy ions provide a nice oportunity to study photon-photon 
collisions [T]. One can expect an enhancement of the cross section for the reactions of this 
type compared to proton-proton or e~^e~ collisions which is caused by a large charges of the 
colliding ions. In this type of reactions virtual (almost real) photons couple to the nucleus as a 
whole. Naively the enhancement of the cross section is proportional to Zl^Z^ which is a huge 
factor. We have discussed recently that the inclusion of realistic charge distributions and realistic 
nucleus charge form factor makes the cross section smaller than the naive predictions. Many 
processes has been discussed in the literature. Recently we have also studied some of them. 

We have discussed production of ii~ pairs [2] heavy-quark heavy-antiquark pairs [3] as well 
as production of two mesons: (I^ pairs [4], 'k^tx~ pairs [5] as well as of DI) meson pairs [6]. 

Here we shall summarize the recent works. 



2. Formalism 

2.1. Equivalent Photon Approximation 

The equivalent photon approximation is a standard semi-classical alternative to the Feynman 
rules for calculating cross sections of electromagnetic interactions [7] . This picture is illustrated 
in Fig. [1] where one can see a fast moving nucleus with the charge Ze. Due to the coherent action 
of all protons in the nucleus, the electromagnetic field surrounding (the dashed lines are lines of 
electric force for a particles in motion) the ions is very strong. This field can be viewed as a cloud 
of virtual photons. In the collision of two ions, these quasireal photons can collide with each 
other and with the other nucleus. The strong electromagnetic field is a source of photons that 
can induce electromagnetic reactions on the second ion. We consider very peripheral collisions 
i.e. we assume that the distance between nuclei is bigger than the sum of radii of the two nuclei. 
Fig. [1] explains also the quantities used in the impact parameter calculation. In the right panel 
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Figure 1. A schematic picture of the colhsion and the quantities used in the impact parameter 
calculation. 

we can see a view in the plane perpendicular to the direction of motion of the two ions. In order 
to calculate the cross section of a process it is convenient to introduce the following kinematic 
variables: 

• X = u}/Ea_, where uj energy of the photon and the energy of the nucleus 

• Ea = 'yAniproton = where Ma is the mass of the nucleus and Ea is the energy of the 
nucleus 

Below we consider a generic reaction A A — )• AAciC2 and later consider different examples 
when ci and C2 are leptons, quarks or mesons. In the equivalent photon approximation the total 
cross section is calculated by the convolution: 

a (AA C1C2AA; saa) = J ^{11 ^ C1C2; = ^/xlX2SAA) dn^7 {xi, X2,h) . (1) 

The luminosity function dn^^ above can be expressed in term of flux factors of photons 
prescribed to each of the nucleus: 

/I 1 
(b) d^biiV {uiM) d'baiV {002, bs) (2) 

The presence of the absorption factor 5^^^ (b) assures that we consider only peripheral 
collisions, when the nuclei do not touch each other i.e. do not undergo nuclear breakup. In 
the first approximation this can be taken into account by the following approximation: 



Si,, (b) = e{h- 2Ra) = 9 (|bi - bsl - 2Ra) . (3) 

In the present case, we concentrate on processes with final nuclei in the ground state. The 
electric field force can be expressed through the charge form factor of the nucleus |2j. 

The total cross section for the AA — )• C1C2AA process can be factorized into an equivalent 
photons spectra and the 77 — )■ C1C2 subprocess cross section as: 

a {AA C1C2AA; saa) = J ^(11^ C1C2; 9 (|bi - bs] - 2Ra) 

N (wi, bi) N {uj2, h2) d^bid^bs^^ . (4) 

Ldl LO2 



We introduce the invariant mass of the 77 system: = y/4ujiuj2- Additionally, we define 
Y = ^ (yci + ^£2) rapidity of the outgoing C1C2 system. Making the following transformations: 

= -^e , W2 = -i^e , (5) 



da;ida;2 — )• dWj^dY, 



d{uJi,UJ2) 
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d{w^^X) 

formula Q can be written in an equivalent way as: 

a {AA CIC2AA- saa) = J ^ hi ^ cics; W^^) 9 (|bi - bal - 2Ra) 
N (ojiM) N {t02, ba) X d^bid^bs^dW^^dF . 
Finally the cross section can be expressed as the five-fold integral: 

a {AA C1C2AA; saa) = j ^{11^ ^l^^l-■,W^^) e (|bi - bal - 2Ra) 
xiV (wi, bi) iV (w2, bs) 27r6„ dhm db^ dby^^dW^^dY , 



(6) 
(7) 



(8) 



(9) 



where hx = ihix + 62x)/2, by = {biy + b2y)/2 and bm = bi — 62 have been introduced. The 
formula above is used to calculate the total cross section for the AA — )• AAciC2 reaction as well 
as distributions in 6 = bm, W-y-y = Mc^c2 and Y[ciC2)- 

Different forms of charge form factors of nucleus were used in the literature. We compare 
the equivalent photon spectra for realistic charge distribution and for the case of monopole form 
factor. A compact formula how the photon flux depends on the charge form factors can be found 
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(10) 



where Ji is the Bessel function of the first kind and q is momentum of the quasireal photon. 
The calculations with the help of realistic form factor are rather laborious, so often a simpler 
formula with monopole form factor is used [8j. 

2.2. Charge form factor of nuclei 

The charge distribution in nuclei is usually obtained from elastic scattering of electrons from 
nuclei P|. The charge distribution obtained from those experiments is often parametrized with 
the help of two-parameter Fermi model [lOj: 



P{r) = Po[^ + exp 



(11) 



where c is the radius of the nucleus, a is the so-called diffiusness parameter of the charge density. 
Fig. [2]shows the charge density normalized to unity. The correct normalization is: po, Am(0) = 

for Au and po,P6(0) = for Pb. 

Mathematically the charge form factor is the Fourier transform of the charge distribution [9j : 



^{q) — / — P (^) sin (Qr) rdr . 
J q 



(12) 



Fig. [3] shows the moduli of the form factor calculated from Eq. ()12p as a function of momentum 
transfer. Here one can see many oscillations characteristic for relatively sharp edge of the 
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Figure 2. The ratio of p the charge distibution to poj the density in the center of nucleus. 
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Figure 3. The nioduh of the charge form factor F^m (q) of the ^^'^Au and ^^^Ph nuclei for realistic 
charge distributions. For comparison we show the monopole form factor for the same nuclei. 

nucleus. We show results for the gold (solid hne) and lead (dashed line) nuclei for realistic charge 
distribution. For comparison we show the monopole form factor often used in the literature. 
The two form factors coincide only in a very limited range of q. 
The monopole form factor [8j: 

= (13) 

leads to a simplification of many formulae for production of pairs of particles via photon-photon 
subprocess in nucleus-nucleus collisions. In our calculation A is adjusted to reproduce root mean 
square radius A = w A with the help of experimental data [lU] . 



2.3. Exclusive production offi~^fj,~ pairs 

Elementary cross section for charged leptons can be calculated within Quantum Electrodynam- 
ics. Several groups have made relevant calculations (see e.g. |1H \T1[ \TE[ [TJ] and references 
therein) . 

Recently we have performed calculation of exclusive production of ^'^ and explored 
potential of RHIC and LHC in this respect. In Ref.[2j we have presented several distributions 
in muon rapidity and transverse momentum for RHIC and LHC experiments, including 
experimental acceptances. We have demonstrated how important is inclusion of realistic form 
factor in order to obtain realistic distributions of muons for RHIC and LHC. Many previous 
calculations in the literature concentrated rather on the total cross section and did not pay 
attention to differential distributions. Hovever, future experiments will measure the cross section 
in very limited part of the phase space. Here we wish to present only some selected examples. 




5 10 15 20 

Pj,, P4, (GeV) 



Q. 
X) 

\ _ 
b 10 



Q. -2 
-O 10 



b 

T3 



10 

o 
I— 

< -4 

10 





1 1 1 1 1 1 1 1 1 1 1 1 1 . 

STAR 






ratio= 


realistic F,^ 






monopole F„ 


■ (b) 





10 15 20 

P3„ p4, (GeV) 



Figure 4. (left) and the ratio (right) for the STAR conditions: 2/3,2/4 G (-1, 1), PuiPit > 1 GeV 
and Wnn = 200 GeV. 

The distribution in the muon transverse momentum for STAR detector is shown in FigUJ 
The STAR rapidity cuts -1 < 2/3,2/4 < 1 are taken here into account. As can be seen from the 
figure, the inclusion of realistic charge distribution is here extremely important. The relative 
effect of damping of the cross section with respect to the results with the monopole charge form 
factor (often used in the literature) is shown in the right panel. At pt = 10 GeV the damping 
factor is as big as 100! Experiments at RHIC have a potential to confirm this prediction. 

The ALICE collaboration can measure only forward muons with psudorapidity 3 < r/ < 4 
and has relatively low cut on muon transverse momentum pt > 2 GeV. In Figl5] (left panel) we 
show invariant mass distribution of dimuons for monopole and realistic form factors including 
the cuts of the ALICE apparatus. The bigger invariant mass, the bigger the difference between 
the two results. The same is true for distributions in muon transverse momenta (see the right 
panel) . 

2.4. Exclusive production of cc and bb 

In Fig l6l7l^nU] we show several photon-photon processes leading to the QQ in the final state. In 
the following we shall discuss them one by one. 

Let us start with the Born direct contribution. The leading-order elementary cross section 
for 77 — )• QQ as a function of W^^ takes a simple form which differs from that for 77 — )• l~^l~ 
by color factors and fractional charges of quarks. 
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Figure 5. 
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(right) for ALICE conditions: 2/3,2/4 = (3,4), P3t,P4t > 2 GeV. 
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Figure 6. Representative diagrams for the Born amphtudes. 
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Figure 7. Representative diagrams for the leading-order QCD corrections. 



In the current calculation we take the following heavy quark masses: rUc = 1.5 GeV, 
nih = 4.75 GeV. It is obvious that the final QQ state cannot be observed experimentally due 
to the quark confinement and rather heavy mesons have to be observed instead. Presence of 
additional few light mesons is rather natural. This forces one to include also more complicated 
final states. 

In contrast to QED production of lepton pairs in photon-photon collisions, in the case of 
QQ production one needs to include also higher-order QCD processes which are known to be 




Figure 8. Representative diagrams for QQqq production. The oval in the figure means a comphcated 
interaction which is described here in the saturation model as explained in the main text. 




Figure 9. Representative diagrams for the single-resolved mechanism. The shaded oval means either 
t- or u- diagrams shown in Fig. [S) 



rather significant. Here we include leading-order corrections only for the direct contribution. 
In a^-order there occur one-gluon bremsstrahlung diagrams (77 — )• QQg) and interferences of 
the Born diagram with self-energy diagrams (in 77 — )• QQ) and vertex-correction diagrams (in 
77 ~^ QQ)- The relevant diagrams are shown in FiglJl We have followed the approach presented 
in Ref. |15j . The QCD corrections can be written as: 



4 2vra2^^^as^(i) 
vr 



77 



(14) 



The function /(^) is calculated using a code provided by the authors of Ref. [TS]. In the present 
analysis the scale of is fixed at pP' = ^tuq. 

We include also the subprocess 77 — )• QQqq, where q (q) are u, d, s, quarks (antiquarks). The 
cross section for this mechanism can be easily calculated in the color dipole framework [16\ I17j. 
In the dipole-dipole approach [17] the total cross section for the 77 — t- QQ production can be 
expressed as: 



= E 

f2^Q 

+ E 



^"^^ (Pi,^i) 



^■^^■^^ (P2, Z2) add {pi,P2,XQf) d^ pidzid^ p2dz2 



{P2,Z2) add {pi, P2, Xfq) d"^ pidzid"^ p2dz2 



(15) 



where ^'^'^ (p, z) are the quark - antiquark wave functions of the photon in the mixed 
representation and add is the dipole-dipole cross section. Eq. (fT5|) is correct at sufficiently high 



energy W^^ ^ 2mQ. At lower energies, the proximity of the kinematical threshold must be 
taken into account. In Ref. [16] a phenomenological saturation-model inspired parametrization 
for the azimuthal angle averaged dipole-dipole cross section has been proposed: 



a 



_a,b 
dd 



1 — exp 



'cflF 
4:Ro (Xab) 



(16) 



Here, the saturation radius is defined as: 



Ro (Xab) 



1 / Xab 

Qo V Xo 



-A/2 



(17) 



and the parameter Xab which controls the energy dependence is given by: 



Xab 



4m^ + 4m^ 



(18) 



The effective radius is parametrized as r^g = {piP2)^ / {P1+P2) ilSj • Some other parametrizations 
of the dipole-dipole cross section were also discussed in the literature. The cross section for the 
77 ~^ QQqQ process here is much bigger than the one corresponding to the tree-level Feynman 
diagram as it effectively resums higher-order QCD contributions. 

As discussed in Ref. jl7j the QQqq component have very small overlap with the single- 
resolved component because of quite different final state, so adding them together does not lead 
to double counting. The cross section for the single-resolved contribution can be written as: 



ai-res {s) = / dxi gi ag^ (s = Xi-s) 



+ 



92 [X2,P- ) Cr^g (S = X2S) 



(19) 



where gi and g2 are gluon distributions in photon 1 or photon 2 and and a-yg are elementary 
cross sections. In our calculation we take the gluon distribution from Ref. |18| . 

Elementary cross sections have been presented and discussed in Ref. [3]. Here we show only 
nuclear cross sections. In Fig. [10] we compare the contributions of the different mechanisms 
as a function of the photon-photon subsystem energy. For the Born case it is identical as a 
distribution in quark-antiquark invariant mass. In the other cases the photon-photon subsystem 
energy is clearly different than the QQ invariant mass. These distributions reflect the energy 
dependence of the elementary cross sections. Please note a sizable contribution of the leading- 
order corrections close to the threshold and at large energies for the cc case. Since in this case 
W^^ > Mqq, it becomes clear that the QQqq contributions must have much steeper dependence 
on the QQ invariant mass than the direct one which means that large QQ invariant masses 
are produced mostly in the direct process. In contrast, small invariant masses (close to the 
threshold) are populated dominantly by the four-quark contribution. Therefore, measuring the 
invariant mass distribution one can disentangle some of the different mechanisms. As far as this 
is clear for the cc it is less transparent and more complicated for the bb production. In the last 
case the experimental decomposition may be in practice not possible. 

In Table 1 we show partial contribution of different subprocesses discussed above. 



2.5. Exclusive production o/vr'^vr^ pairs 

In this subsection we discuss production of only "large" invariant mass vt+tt^ pairs. Brodsky 
and Lepage developed a formalism [19] how to calculate relevant cross section. Typical diagrams 
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Figure 10. The nuclear cross section as a function of photon-photon subsystem energy Wj~^ in EPA. 
The sohd hne denotes the results corresponding to the Born amplitude (cc -left panel and bb -right 
panel). The leading-order QCD corrections are shown by the dash-dotted line. For comparison we show 
the differential distributions in the case when an additional pair of light quarks is produced in the final 
state (dashed lines) and for the single-resolved components (dotted line). 



Table 1. Partial contributions of different mechanisms at ^/sJfN = 5.5 TeV. 
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Figure 11. Typical Feynman diagrams describing the 77 {q<i){<1<i) t^t^ amplitude in the LO 
pQCD. 



of the Brodsky-Lepage formalism are shown in Fig. [TTJ The invariant amplitude for the initial 
helicities of two photons can be written as the following convolution: 



M{Xi,X2)= 1^ dxj^ dy(t>n(x,fil^T^^^^ (^x,y,fi'^^(t>^(y,fiiy 



(20) 



where fix = min {x, 1 — x) \/ s(l — z'^), fiy = min (y, 1 — y) \/s{l — z^); z = cos 9 [19]. We take 
the helicity dependent hard scattering amplitudes from Ref. [20j . These scattering amplitudes 
are different for 7r"'"7r~ and The distribution amplitudes are subjected to the ERBL pQCD 

evolution [211|22]. The scale dependent quark distribution amplitude of the pion can be expanded 
in term of the Gegenbauer polynomials: 



2^/3 



6x {l-x)Y, Cl'^ {2x - 1) an . 



n=0 



(21) 



where /jr is the pion decay constant. 

Different distribution amphtudes have been used in the past. Recently Wu and Huang [23] 
proposed a new distribution amphtude (based on a hght-cone wave function): 



2^/2^3/2/. 



y/x{l-x) (l + Bx dj"^ (2.T - 1)) 



Erf 



^x (1 — x) 



Erf 



2 a; (1 



(22) 



The pion distribution amplitude at the initial scale is controlled by the parameter B. They have found 
that the BABAR data for pion transition form factor at low and high transferred four- momentum squared 
regions can be described by setting B k, 0.6. This pion distribution amplitude is rather similar to the well 
know Chernyak-Zhitnitsky 24 distribution amplitude {((jt^cz ~ 30a;(l — x){2x — 1)^). In the following 
we shah use B = 0.6 and = 0.3 GeV. Then A = 16.62 GeV^i and (3 = 0.745 GeV. 

The total (angle integrated) cross section for the process can be expressed in terms of the amplitude 
of the process discussed above as: 

27r 



4 • 647r2M^2 q 



^ J2 \MiXi,X2)fdz 



(23) 



Ai,A2 



9s \s J 



(24) 



where the factor 4 is due to averaging over initial photon helicities. 

The hand-bag model was proposed as an alternative for the leading term Brodsky-Lepage pQCD 
approach |25j . It is based on the philosophy that the present energies are not sufhcient for the dominance of 
the leading pQCD terms. As in the case of BL pQCD the hand-bag approach applies at large Mandelstam 
variables s ~ — i ^ —u i.e. at large momentum transfers. Diehl, Kroll and Vogt presented a sketchy 
derivation [25j obtaining that the angular dependence of the amplitude is cx 1/sin^ 6. In this approach 
the ratio of the cross section for the tt^tt" process to that for the tt+tt^ process does not depend on 9 
and is ^. The nonperturbative object (s) in the hand-bag amplitude, describing transition from a 
quark pair to a meson pair, cannot be calulated from first principles. In Ref. [25j the form factor was 
parametrized in terms of the valence and non- valence form factors as: 

^TTTT (S) = — a„ y — j 

The a„, n„, and Ug values found from the fit in Ref. 
have averaged these values and used: a„ — 1.375 GeV2, 

In Ref. [5] we have discussed in detail elementary cross sections as a function of photon-photon energy 
and as a function of cos{9). Here we will present only nuclear cross sections calculated within EPA 
discussed in the theoretical section. 

In Fig. [T2]we show distribution in the two-pion invariant mass which by the energy conservation is 
also the photon-photon subsystem energy. For this figure we have taken experimental limitations usually 
used for the tttt production in e^e~ collisions. In the same figure we show our results for the 77 collisions 
extracted from the e+e~ collisions together with the corresponding nuclear cross sections for tt+tt" (left 
panel) and tt^tt" (right panel) production. We show the results for the standard BL pQCD approach 
with and without extra form factor (see [2S]). 

Comparing the elementary and nuclear cross sections we see a large enhancement of the order of lO** 
which is, however, somewhat less than Z\Z\ one could expect from a naive counting. 



slightly depend on energy. For simplicity we 
= 0.4175, as 0.5025 GeV2 and n, = 1.195. 



2. 6. Exclusive production of p p pairs 

At low energies one observes a huge enhancement of the cross section for the elementary process 77 — p^p" 
(see left panel of Fig fT5|) . In the right panel we show predictions of a simple Regge-VDM model with 
parameters adjusted to the world hadronic data. More details about our model can be found in our 
original paper |3| . 

In FigHDwe show distribution in p^p" invariant mass (left panel) and the ratio of the cross section 
for realistic and monopole form factors. 




Figure 12. The nuclear (upper lines) and elementary (lower lines) cross section as a function of photon- 
photon subsystem energy in the b-space EPA within the BL pQCD approach for the elementary 
cross section with Wu-Huang distribution amplitude. The angular ranges in the figure caption correspond 
to experimental cuts. 




Figure 13. Elementary cross section for 77 — > reaction. The fit to the experimental data is 
shown in the left panel and our predictions for the high energy in the right panel. 



2.1. Some comments and outlook 

We have presented some examples of processes that could be soon studied at RHIC or LHC. In all cases 
we have obtained measurable cross sections. We have pointed out that the inclusion of realistic charge 
form factor is necessary to obtain realistic particle distributions. 

Measurements of the processes discussed here are not easy as one has to assure exclusivity of the 
process, i.e., it must be checked that there are no other particles than that measured in central detectors. 
In all cases fissibilty studies, including Monte Carlo simulations, are required. 

In the close future one may expect results for two-pion and single vector mesons production from LHC 
experiments. Exclusive production of one or two pairs of charged leptons should be feasible too. 

Acknowledgments Some of the results presented here were obtained in collaboration with Wolfgang 
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